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. . - We study the general chaotic features of dynamics of the phantom field modelled in terms of a 

single scalar field conformally coupled to gravity. We demonstrate that the dynamics of the FRW 
' model with dark energy in the form of phantom field can be regarded as a scattering process of two 

»' ' types: multiple chaotic and classical non-chaotic. It depends whether the spontaneously symmetry 

(N ■ 

breaking takes place. In the first class of models with the spontaneously symmetry breaking the 

(~| ' dynamics is similar to the Yang- Mills theory. We find the evidence of a fractal structure in the phase 

' space of initial conditions. We observe similarities to the phenomenon of a multiple scattering process 

, around the origin. In turn the class of models without the spontaneously symmetry breaking can 

QQ ■ be described as the classical non-chaotic scattering process and the methods of symbolic dynamic 

' are also used in this case. We show that the phantom cosmology can be treated as a simple 

, model with scattering of trajectories which character depends crucially on a sign of a square of 

^ ■ mass. We demonstrate that there is a possibility of chaotic behavior in the fiat Universe with a 

^\ ' conformally coupled phantom field in the system considered on non-zero energy level. We obtain that 

( \ the acceleration is a generic feature in the considered model without the spontaneously symmetry 

^Nj , breaking. We observe that the effective EOS coefficient oscillates and then approach to = — 1. 

00 ■ 
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I , I. INTRODUCTION 

(D , 

. The observations of distant supernovae [HQ give the evidence that our Universe is undergoing accelerated expansion 
^ ■ in the present epoch. In principle, there are two different approaches to treat this phenomenon. In the first approach 
. it is postulated that there is some unknown exotic matter which violates the strong energy condition p -\- > 0, 
where p is the pressure and p is the energy density of perfect fluid. This form of matter is called dark energy. In the 
^ . past few years different scalar field models like quintessence and more recently the tachyonic scalar field have been 
■ - - ' conjectured for modelling the dark energy in terms of sub-negative pressure p > —p. A scalar field with super-negative 
pressure p < —p called a phantom field can formally be obtained by switching the sign of the kinetic energy in the 
Lagrangian for a standard scalar field. For example in the Friedmann-Robertson- Walker (FRW) model the phantom 
field minimally coupled to a gravity field leads to p + p = —ip^, where p^ — —jip'^ + V{tjj), p^ — — ^V'^ — V^ip), and 
V{ip) is the phantom potential. Such a field was called the phantom field by Caldwell [3] who proposed it as a possible 
explanation of the observed acceleration of the current Universe when flyn,o ^ 0.2. Note that a coupling to gravity in 
the quintessence models was also explored [1]. 

The second approach called the Cardassian expansion scenario has recently been proposed by Freese and Lewis 
as an alternative to dark energy in order to explain the current accelerated expansion of the Universe. In this 
scenario the Universe is flat and matter dominated but the standard FRW dynamics is modified by the presence of 
an additional term p" such that 3H^ = pcs — p + 3Bp^, where H = {d\na)/dt is the Hubble parameter; and a 
is the scale factor. However, let us note that this additional term can be interpreted as a phantom field modelled 
by the equation of state p = p{p) = + 7) — l]/0, where p = Pm.oO'^^^^^'^'' ■ Therefore for dust matter we obtain 
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p = {n ~ l)p, and n < leads to the phantom field. The Cardassian expansion with n < which can be interpreted 
as the phantom fluid effect. 

Phantom fields lead to the super-accelerated expansion of the Universe, i.e. d H/ dt > 0. The simplest models 
describe this field in terms of minimally coupled real scalar field with the negative kinetic energy term — ^V'^ [1, B0]- 
It is interesting that phantom fields are also present in string theories lol. ITol . and arise as a phenomenological 
description of quantum effects of particle production in terms of bulk viscosity [11| . Because the phantom fields violate 
the Lorentz invariance condition most physicists believe that the phantoms open the doors on new physics [l^. The 
investigation of the theoretical possibility to describe dark e nerg y in terms of phantom field was the subject of many 
papers [H [11 [H [H, [13, [H, [S [13, [11, [11 [11 [11 [13 [11, [11, [11 [p . The review on dark energy models 

were presented by Copeland et al. [32|. 

In the paper by Dabrowski et al. [33[ it is considered the quantization of phantom field via the Wheeler-DeWitt 
equation in quantum cosmology. They showed that quantum effects give rise to avoiding the big-rip singularity. Also 
some other basic problems in cosmology like the problem of direction of time can be solved. 

If we confront the phantom field model with the observation of SNIa data we obtain that it is the best candidate 
together with the ACDM model for the substantial form of dark energy [H, [H [s^, [sj ■ 

It is interesting to investigate new dynamical effects created in the presence of phantom fields. The main goal of this 
paper is to model the phantom fields in terms of scalar fields with a potential rather than in terms of the barotropic 
equation of state. For the latter case the dynamics is regular and can be represented on the two-dimensional phase 
space [38| . 

Usually phantom fields are modelled in terms of minimally coupled scalar fields The case of minimally coupled to 
gravity phantom fields was analyzed in the context of existence of periodic solutions [s^. They demonstrated that 
the dynamics is trivial in a sense of nonexistence of periodic solutions. Faraoni used the framework of phase space 
for investigating the dynamics of phantom cosmology, late-time attractors, and their existence for different shapes of 
potential [40] ■ It was showed that dynamics of the flat FRW universe can be reduced to the form of a two-dimensional 
dynamical system on a double sheeted phase space. It is a simple consequence of an algebraic equation for ip which 
can be expressed as a function of the Hub ble p arameter H and tp. In this case there is no place for chaotic behavior 
of trajectories in the phase space {H, ip, tp) [41[. The exit on the inflationary epoch and bounce was showed in the flat 
FRW universe with two interacting phantom scalar flelds [i^ . In the closed FRW cosmological model with minimally 
coupled scalar field there appears transient chaos which has character of the scattering process [11] . The scattering 
process takes place during the bounce — a transition from a contracting to expanding universe. In this context the 
language of symbolic dynamics and topological entropy was used (4^, l45l [46t . 

The conformally coupled scalar field are also of interest. The significance of long- wavelength modes in the WKB 
approximation of a conformally coupled scalar field was analyzed in the inflationary scenario [47] . 

While in this paper we concentrate on the conformal coupled scalar flelds, it is worth to mention the works about 
the case of the generic (non-minimal and non-conformal) coupling between a phantom fleld and gravitation. Faraoni 
\a§^ discusses different motivations for choosing £, ^ (both positive and negative) from the theoretical point of view 
and from observational constraints. The inflation and quintessence with the non-minimal coupling were studied in 
the context of the formulation of some necessary conditions for acceleration of the universe [i^, [lo] (see also [5l|). 
Faraoni also pointed out that the non-minimally coupling term different from the conformal coupling value (C = 1 /6) 
spoils the equivalence principle of general relativity. 

It is interesting that some observational constraints on the non-minimal coupling (possitive or negative) can be 
found from CMB observations. Tsujikawa and Gumjjudpai [52| studied this bounds for a special power-law potential. 
The dynamics with the exponential form of the potential of a scalar fleld was also studied [53] . The two-fleld inflation 
models with negative non-minimal coupling and hybrid inflation are also interesting in the context of the large scale 
curvature perturbation [s^ ] . In turn the constraints on the ratio of the self-coupling and non- minimal coupling constant 
can be obtained during the inflation period [ss!] . 

In this paper we ask what kind of dynamics can be expected from the FRW model with phantom fleld. It is 
well known that the standard FRW model reveals some complex dynamics. The detailed studies gave us a deeper 
understanding of dynamical complexity and chaos in cosmological models and resulted in conclusion that complex 
behavior depends on the choice of a time parameterization or a lapse function in general relativity [sl [H Issj . 
Castagnino et al. [1^ showed that dynamics of the closed FRW models with a conformally coupled massive scalar 
field is not chaotic if considered in the cosmological time. The same model was analyzed in the conformal time by 
Calzetta and Hasi [gI] who presented the existence of chaotic behavior of trajectories in the phase space. Motter and 
Letelier (6l| explained that this contradiction in the results is obtained because the system under consideration is non- 
integrable. Therefore we can speak about complex dynamics in terms of nonintegrability rather than deterministic 
chaos. The significant feature is that nonintegrability is an invariant evidence of dynamical complexity in general 
relativity and cosmology [62|. 

We can find many analogies between models with spontaneously symmetry breaking and the Yang-Mills systems. 
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Problems of chaotic behavior in the Yang-Mills cosmological models have been investigated by Barrow and his col- 
laborators. Chaos from Yang-Mills can only occur in an anisotropic universe but it occurs for some arbitrarily small 
anisotropy (see the study of Bianchi I [63| and other Bianchi types [13] )■ However, when a perfect fluid is added 
things change in an interesting way that contrast with the situation with a magnetic field [65| . 

For the FRW model with phantom it can be shown that there is a monotonous function along its trajectories and 
it is not possible to obtain the Lyapunov exponents or construct the Poincare sections. Therefore it is useful to 
study the nonintegrability of the phantom system and set it in a much stronger form by proving that the system 
does not possess any additional and independent of Hamiltonian first integrals, which are in the form of analytic or 
meromorphic functions. Of course, it is not the evidence of sensitive dependence of solution on a small change of 
initial conditions. However, it is the possible evidence of complexity of dynamical behavior formulated in an invariant 
way [H. 

The notion of deterministic chaos is a controversial issue in the Mixmaster models. The value of the numerically 
computed maximal Lyapunov exponent for those systems depends on the time parameterization used [6^. However, 
the existence of fractal structures in the phase space provides the coordinate independent signal of chaos in cosmology 
as it was shown by Cornish and Levin In particular, they found that the Bianchi IX has a form of chaotic 

scattering. The short scattering periods intermittent integrable motion and evolution of the system is chaotic. It is 
similar to a pin-ball machine. 

The main goal of this paper is to show that dynamics of phantom cosmology can be treated as a scattering process. 
Only if the spontaneously symmetry breaking is admitted this process has chaotic character. Evidences of dynamical 
behavior are studied in tools of symbolic dynamics, fractal dimension, and analytically the Toda-Brumer-Duff test. 



II. HAMILTONIAN DYNAMICS OF PHANTOM COSMOLOGY 



We assume the model with FRW geometry, i.e., the line element has the form 

ds^=-df + a^t)[dx^+f{x){de^+sm^Od^^)], (1) 

where 

{sinx, < X < TT fc = +1 
X, 0<x<(X3 fc-0 (2) 
sinhx, 0<X!iOo fc = — 1 

k — 0, ±1 is the curvature index, < ip < 2tt and < 9 < tt are comoving coordinates, t stands for the cosmological 
time. 

It is also assumed that a source of gravity is the phantom scalar field ip with a generic coupling to gravity. The 
gravitational dynamics is described by the standard Einstein-Hilbert action 

Sg^\ml j dx^g{R-2K), (3) 

where rri^ — (SttG)"^; for simplicity and without lost of generality we assume 47rG/3 = 1. The action for the matter 
source is 

Sph = -l [ dxV^i-g^^'-ij^i:, - + 2Um. (4) 



2 _ 

Let us note that the formal sign of is opposite to that which describes the standard scalar field as a source of 

gravity, where U{'4)) is a scalar field potential. We assume 

U{i]j) = imV^^ + ^AV''* (5) 

and that conformal volume / (Px over the spatial 3-hypersurface is unity. ^ is the coupling constant of the scalar field 
to the Ricci scalar 



i? = 6 - -h - + - (6) 
where a dot means the differentiation with respect to the cosmic time t. 
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If we have the minimally coupled scalar field then ^ = 0. We assume a non-minimal coupling of the scalar field 

The dynamical equation for phantom cosmology in which the phantom field is modelled by the scalar field with 
an opposite sign of the kinetic term in action can be obtained from the variational principle S{Sg + Sph) = 0. After 
dropping the full derivatives with respect to time we obtain the dynamical equation for phantom cosmology from 
variation S{Sg + Sph)/Sg = as well as the dynamical equation for field from variation S{Sg + Sph)/6^ = 

+ SHiJ; - ^ + ^Rip = 0. (7) 
dip 

It can be shown that for any value of ^ the phantom behaves like some perfect fluid with the effective energy and 
the pressure in the form which determines the equation of state factor 

_ -1^^ - ug^) + mam + m-] + e^^(2g + 3^^) ^ p.^ 

Formula ^ differs from its counterpart for the standard scalar field [67] by the presence of a negative sign in front of 
the term ip'^ . 

In equation ([5]) the second derivative (■0^) in the expression for the pressure can be eliminated and then we obtain 
p^ + V-' - iH{i,^) ■ - 2e(6e - - 3^(8e - - C/ W + 2^V'^- (9) 

Of course such perfect fluid which mimics the phantom flcld satisfies the conservation equation 

p^+3H{p^+p^)^0. (10) 

We can see that complexity of a dynamical equation should manifest by complexity of . 

Let us consider the FRW quintessential dynamics with some effective energy density given in equation ([5]) . This 
dynamics can be reduced to the form like of a particle in a one-dimensional potential [68] and the Hamiltonian of the 
system is 

n{a',a)^^+V{a) = 0, V{a)=-p^a^ (11) 

where a prime means the differentiation with respect to the conformal time rj. 

The trajectories of the system lie on the zero energy level for fiat and vacuum models. Note that if we additionally 
postulate the presence of radiation matter for which pr oc then it is equivalent to consider the Hamiltonian on 
the level TC = E = const. Of course the division on kinetic and potential parts has only a conventional character and 
we can always translate the term containing -0^ into a kinetic term. 

The dynamics of the model is governed by the equation of motion ([7]), which is equivalent to the conservation 
condition (jTU]) and the acceleration condition 

^ = -p^(l + 3u.^). (12) 
This equation admits the generalized Friedmann first integral which assumes the following form 

- + _ 3^^2^2 „ 3^iJ(^2y ^ 1^2 _ 1 ^ _ (^3) 

2 2 6 

If we postulate existence of radiation in model then left hand of this equation can be negative. 

Let us consider now both cases of minimally and conformally coupled phantom fields. They can model the 
quintessence matter field in terms of matter satisfying the equation of state p^ = w^p^. 

A. Minimally coupled phantom fields 

For minimally coupled phantom fields (^ ~ 0) the function of energy takes the form 

f - ^ + icb'a - <t>a'f - \m^<^^a^ - - |a^ (14) 

where pcft = ^^V"^ + U{il)), V — —pcsa'^, H — + V{a,ip,tp), (f) = a-ip, (/{tp) = ^m'^ip^ + jXip'^ is assumed. 
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B. Conformally coupled phantom fields 



For conformally coupled phantom fields we put ^ = 1/6 and rescale the field ip (j) = ^a. Then the energy function 
takes the following form for a simple mechanical system with a natural Lagrangian function C = \gai3{(f^)' {(1^)' ~ ^(9) 

£ = \ W + - \rn'<f'a' - - (15) 

In contrast to the FRW model with conformally coupled scalar field the kinetic energy form is positive definite 
like for classical mechanical systems. The general Hamiltonian which represents the special case of a two coupled 
non-harmonic oscillators system is 

n = + Viq) = \{pI+ pD + Ax^ + By^ + Cx^ + Dy^ + Ex^y\ (16) 

where A, B, C, £>, and E are constants. 

In order to study the integrability of dynamical systems we use Painleve's approach. Painleve's analysis gives 
necessary conditions for the integrability of dynamical systems and it is the most popular integrability detector. The 
recapitulation of Painleve's analysis for system (|16p was done by Lakshmanan and Sahadevan ||5!||. They found that 
system (|16l) passes the Painleve test and is integrable in the following four cases 

A = A = 'iK 

A = A ^ K 

A = 8A = 3A 

A = 16A ^ 6A. 



This result is in full agreement with conclusions concerning integrability of the two coupled quartic non-harmonic 
oscillator systems. Therefore, phantom cosmology can be considered as coupled quartic non-harmonic oscillators. 
Of course, Painleve analysis gives necessary conditions for the integrability of dynamical systems. However, there 
exist whole classes of integrable systems which do not possess the Painleve property. It means that the Painleve 
approach gives over-restrictive conditions for integrability. It is obvious that for m? = the FRW phantom cosmology 
is integrable because of the possibility to separate of variables in the potential. Then we have two decoupled quartic 
non-harmonic oscillators. 

The non-integrability of the non-flat FRW model with the scalar field with the potential V{(t>) oc (jP was investigated 
in an analytical way by Ziglin , 70] . For a deeper analysis of integrability in terms of Ziglin and Morales-Ruiz and 
Ramis see [38| . 

It would be useful to compare Hamiltonians for the cosmological model with phantom fields and the standard 
cosmological model with scalar fields. Let us consider that A = A = for simplicity of presentation. Then for both 
models we have Hamiltonians 

Hph = \{-pI-pI) + - a') + \m'a^ct.^ (17) 

and 

Hfrw = \{-pI+pI) + + (18) 

If we add a radiation component to the energy momentum tensor, which energy density scales like pr — Pr.o^^^ then 
both systems should be considered on the constant energy level H = £ = Pr.o or the constant pr.o can be absorbed by 
the new Hamiltonian Ti. = Ti. — p-^^ and TL is considered on the zero energy level. 

Let us concentrate on the flat models to analyze the similarities to the Yang-Mills systems. While the standard 
cosmological model is described by the Hamiltonian 

WFRW = ^(-Pa+4) + ^™W. (19) 

the phantom cosmological model is 



n^^^^i-pl-pD + ^m'^a^cf' (20) 
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The crucial difference between system (flQl) and (|20|) lies in the definiteness of their kinetic energy forms. It is 
indefinite and has the Lorentzian signature for the standard model, and it is positive definite for the phantom model. 
As a consequence we obtain that the configuration space for the standard model is whereas the condition pi-.o + 
^m'^cjPa'^ > determines the domain of the configuration space admissible for motion for the phantom model. Note 
that if < (the model with the spontaneously symmetry breaking) then this domain is bounded by four hyperbolas 
in every quarter of plane {a,(j)). The same situation can be obtained for the model without the spontaneously symmetry 
breaking (m^ > 0) and dark radiation (pr o < 0). 

The flat phantom model with the spontaneously symmetry breaking is well known as the Yang-Mills systems which 
have been analyzed since the pioneering paper by Savvidy |7l| . For this system the Lyapunov exponents were found 
[T^ I and the Poincare sections were obtained [73*] . The spatially flat universes filled with the Yang-Mills fields exhibit 
chaotic oscillations of these fields [63, 64, 65]. 

This system was also investigated by using the Gaussian curvature criterion (t^I . Let us apply this criterion to the 
non-flat phantom model. The potential function takes the following form 

V{a, <j)) = -k{<j)'^ - a^) - ^m^(t)^a^ (21) 

According to this criterion the periodic and quasi-periodic orbits appear in the domains of the configuration space in 
which the Gaussian curvature of a diagram of the potential function is positive. The line of zero curvature separates 
these domains from the instability regions where the curvature is negative. If the total energy of the system E 
increases the system will be in a region of negative curvature for some initial conditions and the motion is chaotic. 

Let us consider the dynamical system in an autonomous form in the form (x*)' = f^{x'^) then we linearized around 
the special solution and we obtain equation 




where builds the Jacobian and 5x^ is the deviation vector connecting the points on two nearby trajectories 
corresponding the same value of parameter t. In our case = a, x"^ = a' , x^ = (j), = (f>' then the local instability 
of nearby trajectories are determined by eigenvalues of the Jacobian matrix, i.e. the effective potential V{a, cj)). The 
eigenvalues are 

Mi,2 = \ [-iH ± V9i/2 + (22) 

7 = 71,2 = \ {-V^^ - Vaa ± ^{V^^ + Vaa? " ^{V^^^Vaa " V^^)) 

Therefore the necessary condition for local instability is that at least one of the eigenvalues is positive [t^. These 
positive values decide about the local instability of trajectories. Note that the negative sign of the Gaussian curvature 
of potential V{a,(j)) is the sufficient condition of local instability [zl, [t^ 

sgn K = sgn(K^^Ka - V^,) < (23) 

We can see that test of negative curvature adopted in our case gives to K < 0. However it is not a sufficient 
condition for the chaos. It should be pointed out that this criterion has a purely local character in contrast to the 
Lyapunov exponent. Moreover the compactness of a region admissible for motion is required for chaos existence. 

It is worthwhile to mention that the Toda criterion is only the measure of the local instability of nearby trajectories 
[t?! . Therefore the presented analysis of the order-chaos transition should be combined with the Poincare sections 
and other deeper indicators of the chaotic behavior. 

III. NUMERICAL INVESTIGATIONS OF CHAOS IN PHANTOM COSMOLOGY 

To make the numerical analysis we need to distinguish the the model without the spontaneously symmetry breaking 
(to^ > 0) and the model with the spontaneously symmetry breaking (m^ < 0). In both cases the system is considered 
on the some energy level. 



FIG. 1: Two samples of trajectories for model class I. Initial conditions ao — 0, — 0.779, (jio ~ 0.75 (left) and 00 = 0.751 
(right), a'o > calculated from the Hamiltonian constraint. 

A. > 

We would like to stress out why we use the nonintegrability criterion instead of standard measure like the Lyapunov 
exponents or Poincare sections. The main reason is that in analogy to Castagnino et al.'s work (soj we can find a 
monotonous function along a trajectory. This excludes the property of the recurrence the trajectories or topological 
transitivity in the standard Wiggins chaos definition. Then F{(j), <f) (or -F(a, a)) is a monotonous function of the time 
parameter along the trajectories and trajectories escape to infinity for arbitrary initial conditions. In this case it is 
impossible to construct the Poincare sections. It is obvious for any curvature and positive values of A, A, and m? . 
The phase space as well as the configuration space is unbounded in this case. In contrast to classical chaotic systems 
there is invariant compact chaotic set for this system. It does not mean the system is non-chaotic as it does not follow 
the Wiggins definition of chaos. The essence of this phenomenon has different nature. The system is oversensitive 
with respect to small change of initial conditions like in chaotic scattering processes. 

The typical case of the model is a class of models without the spontaneously symmetry breaking 

class A: A^B^C = D = 0, S = -1 (or > 0). 

In this case the domain admissible for motion is unbounded. For the classical mechanical systems with chaos there are 
chaotic sets of trajectories on a compact invariant submanifold. However it does not mean that the system with an 
unbounded region admissible for motion is nonchaotic because the essence of chaotic behavior is of different nature. 
The system can be oversensitive with respect to small changes of initial conditions — the key ingredient of chaos. In 
this respect the system is similar to chaotic scattering processes Therefore it seems to be natural to use the 

methods of investigation of scattering processes. 

The class A model does not possess the property of sensitive dependence on initial conditions. Moreover, there are 
no fractal structures of basin boundaries in the phase space. Because the scattering of trajectories in the potential 
well is present, we deal with the nonchaotic scattering process. In more details it will be discussed in Sec. IIVI 

The different chaotic behavior was found in the Einstein- Yang-Mills (EYM) color field in the flat Bianchi I spacetime 
[s^ l- It was considered the model with radiation and isotropic curvature in contrast to the Mixmaster models. It 
was showed that the EYM systems are also an example of chaotic scattering. Because the Lyapunov exponents are 
coordinate dependent they cannot be used for invariant characterization of chaos. For this aim the methods of chaotic 
scattering are more suitable. They are extremely useful in this context because of noncompact phase space (also the 
configuration space) — the major obstacle of the standard analysis of chaos. Let us note that the case of their model 
with flat spacetime is analogous to our second case with the cosmological constant. 

B. < 

The idea of the description of a dark energy field in terms of a Higgs field which creates inertial mass through the 
spontaneously symmetry breaking has been investigated lately (79| . For < we can distinguish two subclasses of 
the model with the spontaneously symmetry breaking for which chaotic behavior can be detected 



class I: A=B = C = D = 0, S > (or < 0), 
class II: A = B = = 0, C > (or A < 0), £; > (or < 0). 
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FIG. 2; Two samples of trajectories for model class II. Initial conditions ao = 0, (/jq = 0.779, (po = 0.75 (left) and ijio — 0.751 
(right), a'o > calculated from the Hamiltonian constraint. 
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FIG. 3: The Poincare section for the model class 1, a — (j> and a' > (upper); a = 0, a' > (bottom). Motion in narrow field 
between two hyperbolas is completely regular. Complexity of behavior comes from motion near the origin of the configuration 
space. 

These are flat models with conformally coupled phantom fields for which trajectories have the property of topological 
transitivity in contrast to the case m'^ > 0. In Fig. [T] and [2] we present sample trajectories for both models and 
evolution for every trajectory is last for the same interval of time. From this simple picture we can initially conclude 
that the systems are sensitive on initial conditions. 

The first class of models is isomorphic with the well known Yang-Mills systems with the potential function V oc x^y^. 
Because the domain admissible for motion is bounded and its boundary has negative curvature, the property of 
recurrence of trajectories is present. We consider the both classes of systems on some distinguished energy level 
H = £ (X pi-fi > 0. The Poincare sections for both classes are represented in Figs. [H IH [5] In Fig. [3] it is shown the 
Poincare section for the fiat cosmological model with vanishing A, A, and < (class I). If we postulate additionally 
the existence of radiation matter in the model than we deal with the system on the constant non-zero energy level. In 
this case we obtain some chaotic distribution of points on the Poincare section (a,(j)'). In Fig. 0] and is presented 
the flat cosmological model with the negative cosmological constant A and < (class II) on the plane {a, a'). In 
this case we can observe islands of stability. However, the trajectories wander to the non-physical region of a < 0. 
This allows many cycles to be considered if we continue the scale factor into negative values. But what it means 
physically is not clear. 

All these figures illustrate what we proved earlier, namely the standard chaos with recurrence of orbits and the 
property of sensitive dependence on initial conditions. 

We cannot perform the analogous analysis for the class of models with > because there is no chaos in 
Wiggins' standard sense. However they possess the property of complex behavior of trajectories similar to the non- 
chaotic scattering process. Note that the evidence of chaos in terms of fractal basins, Cantori or stochastic layers 
requires the recurrence of trajectories. Similarly the Poincare sections can be constructed from many cycles for a 
useful picture to emerge [ssj . 

In the case of < 0, trajectories of the system return to the neighborhood of the origin time after time and we 
have the multiple scattering process on the potential walls (boundaries of the domain admissible for motion) . We can 
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FIG. 4: The Poincare section for the model class II, = and 0' > 0. The empty region at the upper section comes from 
motion of a particle-universe along the (a = 0) axis. 
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FIG. 5: The Poincare section for the model class II, a = and a' > 0. 

control and count how many times the trajectory gets closer to the origin than some fixed distance from the origin. 
With increasing this fix distance the number of controlled scattering events increases and a more complicated fractal 
structure arises. Fig. [HI [7] and [5] illustrate fractal structures of the phase space of initial conditions for trajectories 
reaching some values af and 0/ at some moment of time evolution. The fractal dimension calculated by counting 
cells which contain light and dark areas (i.e. leading to two types of the final outcome) is a/ = = 2 — Dq — 1.676; 
a/ = = 3 — Do = 1.888; a/ = </>/= 4 — Dq = 1.953; a/ = </)/ = 5 — I?o = 1.971. The chosen method of 
counting cell enables us to order fractals with respect to the increasing complexity, i.e. for a nearly regular system 
the fractal dimension is close to one, in turn for a completely chaotic one is approaching two. One can also observe 
that complexity of dynamical behavior of trajectories grows more the longer trajectories come more often in the 
neighborhood of the origin of the coordinate system because the motion between branches of hyperbolas is regular 
(see Fig. ED. 

Fig . [51 [TUl [TT] and 1 1 2 1 present the fractal structure of phase space of initial conditions for the category II model. At 
enough small value of 0/ the system is regular and no chaotic behavior is present (Fig. [9]) . With an increasing value 
of 0/: (j)f = 2 — -Do = 1.812; (j)f = 3 — Dq = 1.967, we observe the transition to chaos which cause to emerge the 
fractal structure of the space of initial conditions (Fig. [TUHT^ . 
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FIG. 6: The fractal structure of the phase space of initial conditions of trajectories chosen at a = and a' > for class I and 
landing at |a/| = 2, 3, 4, 5 (grey) or at \(j)f \ = 2, 3, 4, 5 (white), respectively. Complexity of the fractal structure increases 
with the final state of the trajectories which means that trajectories spend more time in the region near the point (0, 0) in the 
configuration space. The solid areas in the center of figures come from the motion along the (^-axis. 



In Fig.[T3]we plot the phase space of initial conditions chosen at a = and a' > calculated from the Hamiltonian 
which lead to a given property in one cycle of evolution, i.e. to maximal expansion with a positive value of (grey) 
or negative value of cf) (white) - left image, and back to a final singularity again with a positive value of cf) (grey) or a 
negative value of (white) - right image. Therefore if we consider evolution in a physical domain between initial and 
final singularities this system cannot be chaotic. If we prolong its evolution to the non-physical domain a < then 
we obtain chaos. Note that non-intcgrability indicators measure the true intrinsic complexity of the system in both 
cases (38i] . 



IV. PHANTOM COSMOLOGY AS THE SCATTERING PROCESS 



In this section we investigate behavior of a model without the spontaneously symmetry breaking. This is a model 
with ^ = J5 = C = D and S = -1 (m^ 2) on different energy levels £ > 0, £" = and f < 0. 

In the analysis of this case we explore analogy to the classical system in which appears the chaotic scattering. From 
the equation of motion we obtain that for any initial condition trajectories escape to infinity ((/> —^ ±oo a — *■ ±oo as 
time goes to infinity) . For energy £ > Q the configuration space is unbounded and trajectories pass from one quadrant 
to another. On the other hand for £ < motion is restricted to only this quadrant of the configuration space where 
the initial conditions are. 

In Fig.[T3]we present the analysis of dependence of scattering on the energy level. For this aim the initial condition 
in the configuration space are chosen on the line = — a + 10 with a' = (j)' < under the conservation of the 
Hamiltonian constraint. For £ < trajectories go toward the origin of coordinate system and after some time escape 
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FIG. 7: Magnification of the fractal structure of the phase space of initial conditions from Fig. [S] 

to infinity passing through a hne 4> = —a + 10. On the x-axis we mark the initial value of a and y-axis we mark final 
distance from the symmetry point with coordinates (5, 5) at the moment of intersection of line of initial conditions. In 
this figure there is no discontinuities, we do not also observe the fluctuation of distance df. It enable us to conclude 
that we deal with the scattering process although the domain of interaction is not finite |80j . So there is no chaotic 
scattering. The motion of the system is regular and there is no sensitivity of motion with respect to initial conditions. 

For a deeper confirmation of this statement we study numerically a larger number of initial conditions (Fig. llSp . 
In the configuration space the initial conditions are chosen as previously, i.e. (j) = —a + 10 and the initial conditions 
for velocities are parameterized in a natural way by an angle a, i.e. a' = —V{a,4>,£)cosa, <j)' — ~V{a,(j),£)sa\a 
where V{a, 0, £) = \/ E + ni^a^cj)^. In Fig. [15] (a) and (b) it is presented results of our analysis for £ < and £ = 
respectively. The grey area corresponds to initial conditions for which trajectories pass through line </> = — a + 10 with 
(j> > a. Fig. [15] (c) illustrates domains of initial conditions for £ > for which trajectories outcomes reach different 
quadrants : I quadrant - dark grey, II quadrant - grey. III quadrant ~ black, IV quadrant — medium grey. The black 
line presents initial conditions taken from Fig. [T3] We cannot observe chaotic scattering in this case. 

V. ACCELERATION IN PHANTOM COSMOLOGY 

The current Universe is in an accelerating phase of expansion that's why we check whether the scalar fields in the 
FRW cosmology can explain this phenomenon. For this aim let's are consider acceleration equation which can be 
obtained from the Raychaudhury equation 



(X 



(24) 
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FIG. 8: Magnification of tlie fractal structure of the phase space of initial conditions from Fig. [T] 



In the special case of the phantom field minimally coupled to gravity ^ = and without radiation pr = we have 

-= 2(^-2 + [/(V-)). (25) 
a 

Next, for zero energy level and A = from we have that 

H"^ = 2p^ = -V? + 2[/(Vj). (26) 
and finally we receive the acceleration equation which independent of the form of the potential function 

- = 31/^2 + > 0. (27) 
a 

In the case of conformally coupled phantom field ^ = 1/6 we have 

- = -{-^^ - 2[/(V) - ffV' - H{i,^)+i,^{^)) - 2pr (28) 
and after the rescaling time and field variables: dt = ad?7, = aip and inserting the potential function ([5]) we have 

- = -M^-\M'~'iprfl): (29) 

where prime denotes differentiation with respect to the conformal time and dot with respect to the cosmological time. 
We can also express equation of state parameter ([S]) in these new variables 



^(-|0-^-lmW + l^^'^') 
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FIG. 9: The structure of phase space of initial conditions for class II trajectories chosen at = and 0' > and landing at 
(j) = 0.75 (dark grey) and = —0.75 (grey) presented in left figure and (j) = 0.9 (dark grey) and (j) = —0.9 (grey) in right figure. 
The white poles correspond to the centers of stability in Fig. |4l 





FIG. 10: The fractal structure of phase space of initial conditions of trajectories class II landing at <^ = 2 (left) and = 3 
(right). Initial conditions leading to negative values of cj> omitted for clarity. 



We clearly see that this is a good expression only if a 7^ 0. For the physical region a > we can express the effective 
equation of state parameter as 



Weft 



cf> — gm a cf> + ^Prfl 



(31) 



In Fig. [16] we show the evolution of acceleration and evolution of coefficient of the equation of state for model 
class I (flat model with < 0) with respect to the conformal time 77 which is monotonous function of the cosmological 
time in the physical domain. Note that in the case of < the universe is decelerating and stays in the interval 
(—1/3; 1/3), with = —1/3 when 0' = and w^f, = 1/3 when (f) = 0. 

The analogous analysis made for the model class A without the spontaneously symmetry breaking (presented in 
Fig. [TT]) shows that the universe is always accelerating. In this case pressure is always negative. In the case of 
£ > we can have positive value of and this means that density can be negative. Moreover, during evolution 
there are crossings of the line — —1 and then WeB goes to minus one. The amplitude of these oscillations depends 
on pr^o- 
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FIG. 11: Magnification of the fractal structure of the phase space of initial conditions of trajectories class II from Fig. 1101 




FIG. 12: Magnification of the fractal structure of the phase space of initial conditions of trajectories class II from Fig. 1111 

VI. CONCLUSIONS 

In cosmology the Mixmaster models are well known for their chaotic behavior [sHi . In these models chaos is caused 
by anisotropy of space. In this paper we considered the phantom cosmological models which are homogeneous and 
isotropic where the chaos is the consequence of the nonlinearity of the potential function for the scalar field. We 
found chaos is a generic feature of the phantom cosmology with the spontaneously symmetry breaking. It assumes 
the form of the multiple chaotic scattering process. In the case of the absence of symmetry breaking we did not find 
the sensitivity over initial conditions and the system behaves as the standard scattering process. 

There is no difference between the flat universes filled with minimal and conformally coupled scalar fields. As it 
was demonstrated by Faraoni et al. [4l|, [g^ there is no chaos for the former class of models because the solution 
corresponding fc = are restricted to move in some two-dimensional submanifold of the phase space {H, 0, 0). There 
are no enough room for chaotic motion manifestation. The analogous result is valid for the other value of the coupling 
constant ^. Note that we can observe chaotic behavior in some cases of the conformally coupled scalar field in the flat 
universes, when we consider the system on non-zero energy level. Physically it means that the universe is filled with 
radiation matter. 

The reason for which the flat universes with scalar field are non-chaotic lies in lack of invariant measure of chaotic 
behavior in the general relativity [sil Is^ . In general the standard chaos indicators like Lyapunov exponents or 
Kolmogorov-Sinai entropy depends on time parameterization. In general relativity there is freedom of choice of 
the lapse function which defines reparameterization of time. The Hamiltonian is given modulo the lapse function. 
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FIG. 13: The structure of phase space of initial conditions for model class I chosen at a = with a' > for trajectories leading 
to a' = (or equivalent H — 0) with 4> > (grey) and < (white) - the left image, and initial conditions of trajectories 
leading to a = with (f> > (grey) and < (white) - the right image. There is no chaos in one cycle of evolution. 

Therefore if we consider the Hamiltonian FRW system with scalar field on zero energy level it is possible that in the 
Hamiltonian constraint a i a appear only in the combination H = a/a. Hence the motion of the system takes place in 
the 3-dimcnsional phase space h(H, (j),4>) = on some 2-dimensional submanifold. For any value of coupling constant 
and general class of potentials there is no place for chaos (even the models with the cosmological constant). 

The new picture appears if the system is described on a non-zero energy level. Then the motion of the system 
is in the 3-dimensional submanifold of the 4-dimensional phase space. But in general case there is no possibility of 
choosing the lapse function in such a way that a i d contribute in the Hamiltonian constraint by the Hubble function. 

Among the conformally coupled models there arc in principal two types of behavior. First, chaotic scattering 
process takes place for case < (the spontaneously symmetry breaking). In this case trajectories possess the 
property of topological transitivity which guarantees their recurrence in the phase space and the scattering process 
takes place around the origin. It is similar to the chaos appeared in the Yang-Mills theory. Exploring this analogy 
to this system we calculate a Gaussian curvature of the potential function which measure local instability of nearby 
trajectories. We proof that this curvature is negative. Because the configuration space is bounded in this case the 
mixing of trajectories is observed. 

Second, in the opposite case of m? > we characterize dynamics in the system in terms of symbolic dynamics. 
No fractal structure in the space of initial conditions was found and conclude that the scattering process has not the 
chaotic character. 

We also found that the universe filled with conformally coupled scalar field is accelerating without the positive 
cosmological constant. In principle, there are two sources of this acceleration. The acc:c;lc!ration is driven by negative 
energy density and coefficient of equation of state is greater than —1/3. Alternatively, energy density is positive 
and the coeflacient of equation of state can be smaller than —1/3. The latter there is crossing of = —1. Of course, 
for the explanation of SNIa data apart from the acceleration of the Universe itself, the rate of acceleration is also 
crucial. 

The system with phantom scalar field coupled conformally to gravity in the FRW universe can be treated as 

scattering process with chaotic and nonchaotic character. The generic feature of this class of systems with the 
spontaneously symmetry breaking is the chaotic scattering of trajectories. 
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FIG. 16: The evolution of acceleration from Eq. (|29p with A = (left panel) and equation of state parameter from Eq. (I30p 
(right panel) for model class I with respect to conformal time for short period of time for trajectory from left panel of Fig. [T] 
Shaded area denote unphysical regions of scale factor a < 0. In physical regions acceleration is always negative a < model 
decelerates and equation of state parameter never crosses ui^ = — 1 barrier. 
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FIG. 17: The evolution of acceleration from Eq. p9p (left panel) and the equation of the state parameter WcH from Eq. (I3ip 
(right panel) for model class A without spontaneously symmetry breaking for a sample trajectory starting from the same initial 
conditions in the configuration space (ao = 4.75, 0o = —0,0 + 10 and Oq = 0o < calculated from the Hamiltonian constrain) 
but for different energy levels: £ < (top), £ = (middle) and £ > (bottom). For every sample we have acceleration but 
the equation of the state parameter oscillates and aproaches w^b = — 1. 
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